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In the limit of high vibrational and electronic principal 
quantum numbers a semi c I ass i ca 1 model for au to I on i za 1 1 on of 
diatomic molecules has been constructed. The molecular vi- 
brations are treated as classical oscillators, whose parameters 
are chosen from a knowledge, which can be taken either from 
theory or experiment, of the positions of the quantal vibra- 
tional levels and the position and depth of the minimum of 
the adiabatic potential curve of the residual mojecular ion. 

The theory is designed for the case'where vibrational quantum 
numbers are greater than about 10. Electronic quantum numbers 
are restricted by the requirement that they be high enough 
so that the excited eJectron-core electron Interaction can 
be represented by the monopole term In the region where the 
distance of the excited electron from either nucleus is always 
greater than the distance of the core electron and the excited 
electron-nuclear Interactions can be represented by the monopole 
terms In all regions. An advantage of the theory over the per- 
turbed stationary state theory Is that its validity extends into 
the region of very high electronic quantum numbers (n»100) where 
the electron and nuclear velocities are comparable, and the 
Bo rn-Oppenhe I me r theory is not valid. Nume r I ca I ' es 1 1 ma tes for 
the autoloni zat ion rates are presented for several sample cases 
for vibrations In the neighborhood of n'»10 and excited electrons 
of n«10 to n«20, and for zero energy ejected electrons. 
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I n t roduc t 1 on 


A molecule can autoionlzc when one of Its electrons Is so 
highly excited that a transition from an excited vibrational or 
rotational level can supply enough energy to Ionize the electron. 
The electron is almost decoupled from Its parent molecular Ion, 
and the process Is almost the process of energy transfer between 
colliding systems. Typically a n'*5 to n'«4 vibrational transi- 
tion supplies enough energy to lonlzd'an n«9 electron. The pro- 
cess is observed experimentally in the neighborhood of threshbld 
In the photoionization of 1^, in which there is photoexcitation 
of electrons to around n« 1 0 and vibrations to n ' *5 and then the 
production of photoelectric current by au to I on I za t i on . In this 
region of energy an n«10 electron has an order of magnitude 
h I ghe r ve 1 oc i ty than the relative nuclear velocity of a n'«10 
vibration (estimated from the exact vibrational energies). 

The mechanism responsible for the transition Is taken by Ntelsen 
and Berry (hereafter referred to as NB), In an application of 
the perturbed stationary state theory for heavy particle colli- 
sions^, to be the nonadiabatic (vibronlc) corrections on the 
adiabatic electronic states for the nuclear motion; thus the 
transition amplitude Is proportional to vj£iJij(r,R) and (r ,K) • 

where ij/. and Xj are the Initial electronic and vibra- 
tional states and r and ft are the position vectors of the elec- 
tron In the molecule fixed frame and the re 1 a 1 1 vejnuc 1 e I respec- 
tively, Integrated over the product with the final electronic 
and vibrational states. The first derivative term Is dominant, 
and auto I on I za 1 1 on rates calculated by this theory are In 
qualitative agreement with the experimental estimates of Chupka 
and Be rkowi tz . * 
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The purpose of this paper Is to formulate a theory In 
which a u t o I on i z a t I on rates can be predicted without an explicit 
knowledge of the vibrational states Involved, which In Its turn 
requires a knowledge of accurate adiabatic potential curves and 
surfaces, available for only the smallest molecules. Provided 
the vibrations are excited enough (about n'-10) so that they 
can be represented by classical oscillators, the Information 
required about the vibrational levels involved !n the Internal 
conversion process can instead be taken from the extensive 
infrared spectroscopic measurements of vibrational quanta and 
experimental determinations of molecular b I hd i ng- energ Ics . This 
semi emp i r I ca 1 aspect should be especially useful for polyatomic 
molecules (not explicitly considered here). In addition the 
formulation Is not made within the framework of the Born-Oppenhe I mer 
theory and thus avoids many of the numerical complications of this 
theory and more importantly extends into the region of comparable 
e 1 ec t ron - n u c 1 e a r velocities in which the BO theory is not valid. 

In the energy range considered by NB It Is reasonable that non- 
adiabatic effects could be irfvoked to give a reliable first-order 
theory. However In the region In which SOcnccOthe electron and 
nuclear velocities are comparable so that vibronic couplings are 
no longer perturbations on the Born-Oppenhe i me r system. A similar 
situation Is encountered in H atom-H atom collisions above 1 keV 
incident projectile energy In which the ground state-ground state 
He I t 1 e r-London energy curve departs from its zero velocity limit 
below IkeV, distorts beyond recognition from 1 to 100 keV, and 
merges Into a coulomb curve when the electron exchange interactions 
become negligible above 100 keV 

Au to i on i za 1 1 on from the higher levels (electronic and vibra- 
tional) may be Important in HI regions in interstellar space in 
which competing processes for electron supply are ionization of 
n-100 hydrogen atoms by collision with ground state hydrogen atoms 
and by background radiation from H2 regions. There has been recent 
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as t rophys I ca 1 Interest** In "anomalous" recombination radiation 
from H2 regions as a result of observation through HI regions 
In which such processes are occurring. Molecular autoionization 
has not yet been proposed, although undissociated H 2 molecules 
are known to be present in dark dust clouds in HI regions^ 



Theory 


Consider the molecule in an highly excited electronic 
and vibrational state. We choose to treat the vibrational levels 
of the H ^ residual ion for n '29 as classical oscillators. This 
Is a good approximation for the higher states since the classical 
distribution function can represent an average of the oscillatory 
quantum distribution over the stretch of the molecule^. Hence 
a knowledge from experiment or theor'y' of the positions of the 
quantum states and the depth, D , of the minimum of the adiabatic 
potential curve of the residual ion at R g can be used to choose 
the parameters of the classical oscillator,, whose position and 
velocity are given by, 


R|(t) - R^ + DjSinwt 

v j ( t ) ■ D j ojeostut ( I ) 


l 2 2 

where ^“(k/y) 1 where k is the stretch constant equal to (3 V(R)/9R ) 

and V ( R ) is the potential cu,/ve. The stretch constant is taken 

2 f 9 

from the Morse model to be 2D fi a where a is a fitting parameter . 

U Is the nuclear reduced mass. D. is the amplitude of the i 

1 .2 

oscillator and is chosen by fitting the classical energies, ikDj, 
to the quantum eigenvalues^ where the bottom of the well 9 D 9 Is 
chosen as the zero point. We use the Wentzel^ correspondence 
relation which gives a classical velocity which is an average 
of the Initial and final velocities resulting froma quantal 
absorption of energy, 
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v 


(v ( v f (V| + v f )/2) 


h 


( 2 ) 


This "symme t r I zat I on" seems especially justifiable as we approach 
the top of the potential well where the quantum states are almost 
degenerate. We do not symmetrize the trajectories, R ( , but choose 
the Initial trajectory. In a more careful treatment the symme- 
trized trajectories could be obtained from an Integration over 
the time of the correspondence given'Tn (2). 

The excited electron Is In a state of high principal quantum 

number so that Its average distance from the H* residual core 

2 + 

Is large (of order n ). The existence of in a highly excited 
vibrational state means that the quantum mechanical average 
relative position of the nuclei can change over a very wide 
range of Internuclear distances as determined by the modulus 
| Xi (fi) I of the vibrational wave function in distinction to 
the range for a low vibrational state whose modulus Is peaked 
near R fi . Classically this means that there Is a ( large vibra- 
tional amplitude, D., and that in the time, (t 2 -t^)/2, of a 
half period in which the molecular Ion goes from Its minimum 
to maximum stretch the sphere swept out by the nuclei In the 
space occupied by the excited electron has a large radius; 
hence the average electron distance from the core Is effectively 
decreased (the electron, originally outside the core at 
has a greater chance of being inside the core where the Ionizing 
interaction occurs over the period T«t 2 “tj). At t^ we take the 
electron to be outside this sphere (outside the core, for which 
r>R) , so that it moves in the potential, 


1/r - R 2 /r^P 2 (r*K) - (higher order poles) 


( 3 ) 


whers the monopole term of the expansion In the region r>r' 

( r • Is the core electron coordinate) of the excited electron* 

core electron potential r^j has been cancelled (complete 

screening) and the higher order poles dropped. The quadrupole 

term In (3) was first proposed- by Berry'' to be the dominant 
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Ionizing Interaction, but calculations by Berry , Bardsley , 
and Russek e£ a_^' ^ showed that It contributed only a few per- 
cent to the rates as measured by Chupka and Berkowltz'. Physi- 
cally this Is clear because t'e electron Is always outside the 
molecular core (r>R), and vlbronlc effects are of higher order. 
We will there neglect the quadrupole term. Thus at t*t^ the 
electron moves In a coulomb potential, and Its states a.e hydro-, 
genic. In the Interval tp , the period of the vibration, 
the nuclei move from their minimum to maximum stretch and back 
to their minimum. The ionizing interaction occurs during this 
Interval. Also, during this interval It Is reasonable to treat 
the electron as a molecular electron, whose molecular orbital 
tj/j (In the simplest "Linear Combination of Atomic Orbitals" 
approximation) is given by, 


*1 - a 1 (t) (T A )e 


I V»"f ft \ ~ 1 V."?v - I ( 

1 *mm (? B )e )e 


£ „*v 2 /2) 


(M 


where the states are bound hydrogenlc states of energy on 

protons A and B, and the exponential factors result from the 

translational component of motion of the electron travelling 

with the classical protons moving In the directions of positive 

1 k 

and negative velocities respectively The upper and lower 

signs are taken for the gerade and ungerade molecular symmetries. 
The final ipolecutar orbital Is* 


*, • 




i v • r 


^kim ^ r B^ 


^lv’r^-1 (k 2 /2*v 2 /2)t 


(5) 


where the states are the hydrogenlc continuum states of energy 
k /2 on protons A and B. The. total wave function is the sum of 
(M and (5) (in the two-state approximation 1 ^). Substitution 
into the time-dependent Schrodinger equation and projection 
from the left with the molecular orbitals yields the usual set 
of coupled first-order d I f f e ren 1 1 a 1 * fequa 1 1 ons In the time for 
the set of a ( t) ' s , 

S a - < Iti) “ 1 V a (6) 


where $. and V. are the overlap and potential matrices respectively. 
Ourlng the time of ionization, the potential is Just — 1 / K ( t ) + 1/r 
in the region r<R(t), r>r' and zero In the region r>R(t), r>r' 
(zero In the region r>R(t), r>r' is consistent with the neglect 
of the quadrupole term in (3)), where we have kept only the mono- 
pole termf(as in (3)) of |? - 1 and | "r - r'| ' in the region 

r<R(t), r>r', and the monopoles have cancelled in the region 
r > R ( t ) , r > r ' . Note that the excited electron is always taken 
to be in the region r>r'. This bears closer examination since 
r<R(t) implies that the excited e 1 ec t rorjcar. move in the space 
of the core electron during the collision. Hence the Integration 
over r extends only to R(t). We solve (6) by neglecting back- 

i c 

coupling and distortion , thereby reducing the set to the Born 
amp 1 I tude , 


a(T,k) - (I*)" 1 \dt e l(k /2 ‘S»* *(<$ k (^) | 1/r A - l/R ( t ) | * (^) > 



1 < *iam ( V l ( 1 /r A‘ ,/R(t))e 


(7) 


2 I v 


I ^n£,m ^ r > ^ 


where we have reglected the continuum-continuum overlap. The 
momentum factors vanish by orthogonality except for the Inter- 
ference term which results because the exchange Integral contains 
the atomic orbital components of opposite velocity direction 
of the molecular orbital (cross term). Since the time of a 
classical period is 2 tt/u , and the mofiplus of a (T , k) Is the 
probability, which cannot be greater than unity, the seml- 
elasslcal au to I on I za 1 1 on rate has an upper bound of (2it/w) 1 
- 0.65865 x 10 8 sec" 1 . 
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Numerical Results and Discussion 


To obtain numerical estimates of a(T,k) we make the following 
set of approximations. Since the electronic states are continuum 
and highly excited bound states, r^r p to a good approximation, 
i.e. the two-center Integration can be performed about a single 
center. We expand the continuum function In partial waves and 
keep only the K - 1 wave (J(-1 bound states are the Initial states 
of Interest in the pho to i on i za 1 1 on experiment'; henceA?-0 only 
are considered). The maximum velocity, k , of the ejected 

ma X 

electron is small for the small energy transfers Involved, so 
that we evaluate the rates for zero energy electrons only; the 
radial coulomb wave function normalized on the energy scale 
(coulomb functions normalized on the k scale multiplied by 

■.**> .,' 6 . 


♦ot- < 2 ' r > lj 


2 1 + 1 


((Sr) 4 ) 


(3) 


V 

% 


Using (8) the rate per unit velocity Is given by |a(T,0)| /T 
at zero velocity. Assuming the probability to be nearly cons- 
tant over the range Oikik , the total rate Is estimated by 

rn d ^ 

taking the product | a ( T , 0 ) | k max /T* Numerical difficulties 
In evaluating the generalized Laguerre polynomials limit the 
calculation to about n-^20; continuum waves of very *ow k are 
likewise subject to the same numerical difficulties (convergence 
of the confluent hype rgeome t r i c function for large n«-'/k). 

Table 1 gives sample values for the rates for n*"9, 10, 

11, and 12 vibrations (all possible An') and 10^n^«20. Some 
of these rates violate the upper bound given in the previous 
section because the probabilities are greater than unity. This 
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error is inherent in the Born app rox l ma t i on , which does not 

conserve probability. It Is strongly suspected that the larger 

than unity probabilities are due mainly to the overestimation 

of the two-center integral by the one-center approximation and 

(probably a more Important source of error) the treatment of 

the probability as a constant over the interval 0<k£k , 

rUtcs , 

especially since only the An‘>1^show the violation of the 

bound and k increases with An'. This increase of rate with 
max 

increase of An 1 i s also not in agreement with the NB rates, which 
decrease with An'. This incorrect beh'avior of the sem 1 c 1 as s I ca 1 
rates seems to be due (with smaller contributions from the two 
sources of error above) to the dependence only on the Initial 
trajectory; with increase In An' the final trajectory (amplitude) 
has less and less stretch while the Initial trajectory remains 
the same- Sy mme t r i za t i on of trajectories would therefore be 
indicated. The decrease with increase of n is in agreement with NB. 

g 

The upper bound of nearly 10 means that the sem I c l a s s I ca 1 

1 2 

rates can never be as fast as the quantum rates of around 10 
for An * ■ 1 for vibrations O^n'^5 and n^lO of Berry and Nielsen 20 . 
Their values forA^'>l are as slow as 10^. Clearly the very fast 
rates are the result purely of the quantum effect of initial 
and final vibrational states whose overlap is strongly peaked 
in a very small region of i nternucle r distance. This effect 
could be mimicked classically If we allowed 0 <gjT^2tt; that Is 
if we restricted the classical time to some very small region 
instead of allowing the ionization to occur over the entire period. 
Since there is no unique prescription for choosing such shorter 
times, the classical description of the nuclear motion must fall 
for these cases. 

There is an intermediate region (5*0^12) Into which the 
calculation of NB does not extend, so that, a direct comparison 
between quantum and semi c 1 ass I ca 1 rates cannot be made. It would 
be useful if calculations based on the NB theory were made In this 
region. 
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